Emergence of spatial spin-wave correlations in a cold atomic gas 
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Rydberg spin waves are optically excited in a quasi-one-dimensional atomic sample of Rb atoms. 
Pair-wise spin-wave correlations are observed by a spatially selective transfer of the quantum state 
onto a light field and photoelectric correlation measurements of the light. The correlations are inter- 
preted in terms of the dephasing of multiply-excited spin waves by long-range Rydberg interactions. 



The realization of model systems of many particles 
with precisely controlled and widely tunable interactions 
is of considerable interest to improving our understanding 
of many-body quantum physics. Cold atomic gases have 
become a fruitful platform for such studies as they facili- 
tate experiments under well- understood conditions [l|, Q ■ 
Collisional interactions based on low-energy s-wave scat- 
tering have been used in such studies, with their strength 
tuned by Feschbach resonances 0, 0j- Atoms must be 
cooled to ultra-low temperatures for these short-range 
interactions to dominate the dynamics. 

In contrast, Rydberg- level atoms possess long-range 
interactions which scale as 1/i? 6 in the van der Waals 
regime and 1/i? 3 in the dipole-dipole regime. The 
strength of these interactions for values of the princi- 
ple quantum number n > 100 can be comparable to the 
strength of the Coulomb interaction between ions. This 
leads to Rydberg excitation blockade, in which a sin- 
gle excited atom can suppress the excitation of multiple 
nearby atoms [Hl-[l3|. The multi-atom blockade mecha- 
nism has attracted considerable interest for studies at the 
interface of atomic and condensed matter physics HI- 21 1 . 

The interactions between atoms are switched on and 
off using laser fields tuned to the frequency between the 
Rydberg level and the ground level. In an ensemble of 
atoms, such coherent driving generates collective Ryd- 
berg spin waves. Spin waves can be mapped onto light 
using quantum state transfer techniques. Subsequent 
photoelectric detection of the light permits one to extract 
their properties [j^ - Erjj l . Crystals composed of many such 
spin waves can be generated dynamically by an adiabatic 
sweep of laser detuning through the so-called Rydberg 
staircase, with the number of excited atoms determined 
by the sweep parameters 



Recently, Ates and 



27 



Lesanovsky have shown that an atomic gas undergoing 
continuous-wave, on-resonance Rydberg excitation can 
develop long-range correlations arising from the interplay 
of atomic interactions and dynamic spin- wave generation 
211. 



Here we report observations of interaction-induced lo- 
calization of Rydberg spin waves optically-excited in an 
ultra-cold atomic gas. The experimental apparatus and 
methods have been described in Ref. [29j , where we pro- 
duced single Rydberg spin- wave excitations generated in 
a ~ 15 /im-long ensemble. In this work we employ sim- 



ilar techniques to study pair-wise spin- wave correlations 
in a longer, up to ~ 100 /zm, ensemble. 

A sample of 87 Rb atoms of peak density po ~ 10 12 
cm~ 3 is prepared in an optical lattice. The lattice is 
shut off, and an ensemble is driven in two-photon res- 
onance between the ground \g) = |5S]y 2 ) and a Ryd- 
berg |r) = |102si/ 2 ) levels with a laser excitation pulse 
of duration r = 200 ns, Fig. 1. The transverse size 
(Gaussian waists w x w y ~ 6 /im) of the Rydberg ex- 
citation region is determined by the overlap of the two- 
photon excitation laser fields f2i at 795 nm and fi 2 at 
475 nm. The longitudinal size (waist w z variable be- 
tween 10 and 60 ^tm) is determined by the sample size 
along z. The number of atoms N involved in the ex- 
citation is therefore (1 — 5) x 10 3 . The values of £li 
and f2 2 are chosen such that the multi-atom pulse ex- 
citation area 9 = ^/NQt ~ 1, where fl = £!ifi 2 /A is the 
single-atom, two-photon Rabi frequency, and A is the 
laser detuning from the intermediate level |e) = |5pi/ 2 ). 
Since w x and w v are smaller than the effective atomic 
interaction range, quasi-lD spin waves can be realized 
for w z > w x ,w v . In the limit of instantaneous ex- 
citation (r — > 0), the spin- wave can be described by 

a coherent state l^o) = J2a=o( Ca / y ^-)^i^ 1^)' wnere 
the creation operator for the collective excitations is 
S l = (l/v^V) E M e-^-^a]? and &%> is the raising op- 
erator for the atom [i = [1, N], and \G) = YliLi i s the 
collective ground state [22j. Here k = ki + k 2 , and ki j2 
are wave- vectors of the excitation laser fields Oi j2 . In our 
measurement the intensity of the excitation is chosen so 
as to create on average approximately one excitation in 
the sample 22, 2{|. Correspondingly a Poissonian distri- 

/ v/eal, is used for the initial 



bution of unit mean, c a = 1/ V 
spin-wave excitations in the numerical simulations. For 
finite values of r, Rydberg-level interactions suppress the 
excitation of nearby atoms, resulting in a reduction of 
amplitudes c a for a > 2. 

The collective multi-atom state evolves during the en- 
suing storage period of duration T s . For T s shorter than 
the motional dephasing time the spin- wave dynamics are 
dominated by the interactions of Rydberg atoms. As a 
result, the coherently excited sample evolves into a set of 
localized spin waves, each one containing no more than a 
single excitation. This may be viewed as spin-wave anti- 
bunching, in which the probability of having excitations 
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FIG. 1: (color online), (a) Experimental geometry: an atomic sample of temperature T ~ 10 jiK is produced in an optical 
lattice formed by a retro-reflected 782 nm laser beam propagating along the y-axis. The quantization axis is defined by a 
4.3 G bias magnetic field along the a;-axis. Localized spin waves are formed in the sample as a result of Rydberg excitation 
blockade and spin-wave dephasing. The spatial structure is revealed by scanning the focused waist of the probe beam along the 
z— axis of the two-photon excitation region. Either a single (orange) or double (blue) excitation is formed in an approximately 
Gaussian density profile p of the atomic sample. The probe beam destroys the phase-matched spin waves in a spatially-selective 
manner, affecting the statistics of the retrieved field. The latter is measured by a Hanbury Brown- Twiss type setup consisting 
of a beamsplitter and two single-photon detectors Di and D2. (b) Experimental sequence of spin- wave creation, probing, and 
readout. 
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FIG. 2: (color online). Second-order intensity correlation 
function gr' 2 '(0) and the average number of spin-wave exci- 
tations rj as a function of the sample length. The measured 
values ofg (2) (0) are shown as squares, whereas the circles and 
diamonds show calculated g^ 2 \0) and rj, respectively. In the 
numerical simulations the interval between the excitation and 
the retrieval T s — 0.8 fj,s and the transverse sample waists 
w' x y =6.4 /jm are chosen to approximate the conditions of 
the experiment. The dashed line at 1/e shows rj when all spin 
waves with more than one excitation have dephased [22| . 



in two different regions of the sample is suppressed for 
nearby regions. 

The spin wave is converted into a retrieved light field 
by a 1 /is long read-out field ^3 at 475 nm, in resonance 
with the 1 102s) ■<-» |5pi/2) transition. The retrieved field 
is coupled into a single mode fiber followed by a beam 



splitter and a pair of single-photon detectors Di and D2. 
The photoelectric events on the detectors are recorded 
within gated time intervals determined by the length of 
the retrieved pulse. Photoelectric detection probabilities 
for Di and D2 are calculated as pi.2 = Ni^/Nq, where 
Ni,2 are numbers of recorded events and No is the num- 
ber of experimental trials. Photoelectric detection proba- 
bility for double coincidences P12 is calculated as N12/N0, 
where N12 is a total number of simultaneous clicks on 
both detectors. The second order intensity correlation 
function at zero delay is given by g 2 (0) = Pn/ijPiPi)- 
Using weak classical light pulses within the experimental 
sequence we have measured g 2 (0) = 0.99(2), consistent 
with the expected unity value. 

For a sufficiently small sample (w z < 20 /im in 
Fig. 2), the combination of the excitation blockade and 
interaction-induced dephasing can completely suppress 
amplitudes for two or more excitations in the retriev- 
able spin wave, realizing a source of single photons for 
which g^(0) — > (29[. The atoms perform the role of a 
quantum filter: after the (weak) coherent state of light 
is mapped onto atoms, Rydberg interactions dephase the 
components of the spin wave that contain more than one 
excitation. The phase-matching-based read-out protocol 
therefore generates a retrieved field containing the vac- 
uum and the single-photon components only, for which 
the average number of spin- wave excitations rj = 1/e [22| . 

For a longer ensemble, the l/R 6 scaling of the interac- 
tion strength in the van der Waals regime can result in 
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FIG. 3: (color online). g^ 2 '(0] (top panels) and P (bottom 
panels). Left panels are for w z = 42 fim, right panels for 
w z = 64 fim, with these values extracted by the fluorescence 
imaging of the atomic cloud. The solid curves are based on 
computer simulations with w' z = 42 (jm for left panels, and 
w' z — 64 /mi for right panels. The shaded regions are for a 
range of w' z of 38 to 46 /jm, and 58 to 70 /mi for the left and 
right panels, respectively, while the results are only weakly 
sensitive to the variation of w' x and w' y . 

well-defined "interaction volumes" within the sample. In 
this case, the emergence of a spatial pattern of retriev- 
able spin-waves is expected. In Fig. 2 we display the 
measured values of </ 2 )(0) as a function of the ensemble 
length, together with the results of the numerical simula- 
tions of g( 2 >(0) and of the average number of retrievable 
spin-wave excitations r/ = (Sk Sk ) for T s — 0.8 /is. The 
latter should be compared to its initial value of 1 for the 
Poissonian distribution of unit mean. Both gW(0) and 
X] rise with sample length due to the increasing survival 
probability of double excitations. 

To look into the spatial structure of the spin wave, 
we illuminate it with a near-resonant 795 nm probe laser 
field fl p for 150 ns during the 0.5 /is time interval between 
the Rydberg excitation and readout stages, Fig. lb. The 
beam, focused down to a rj 10 /mi waist (1/e 2 intensity) 
at the position of the sample, produces a local dephasing 
of the spin wave via the associated inhomogeneous light 
shifts and spontaneous emission. The retrieved signal in 
this case corresponds to the spin- wave content of the part 
of the sample unaffected by the probe beam. By scanning 
the position of the probe beam z', it is, therefore, possible 
to reveal spatial spin- wave correlations. 

The measured 3^(0) and the probability of photoelec- 
tric detection P = pi + p 2 are shown in Fig. [3] When 
the beam is sufficiently far from the sample so as not to 
affect it (| z' \ > 0.2 mm), the measured values are consis- 
tent with those for no probe beam (the blue diamond in 
each panel) . As the beam position approaches the sample 
notable decreases in g( 2 )(0) and efficiency are observed. 



P reaches its minimum value with the beam illuminat- 
ing the center of the ensemble, while g^ 2 \0) has a peak 
whose prominence depends on the length of the cloud. 

Qualitatively these features can be understood as fol- 
lows: the single excitation is peaked at the center of the 
ensemble because of the higher atom density, while the 
double excitations are localized at the ends of the ensem- 
ble because of the interaction-induced dephasing. The 
two-photon detection events correspond to two spatially- 
localized singly-excited spin waves on the two ends of the 
sample. When the probe beam is focused on one of the 
ends, the corresponding spin wave is destroyed, thus sub- 
stantially reducing the probability of a two-photon event 
P12 ■ The reduction of the numerator results in a low value 
of <? 2 (0). When the probe beam illuminates the center 
of the ensemble, pi2 is not significantly affected because 
the corresponding excitations are localized largely at the 
sample ends. The probabilities p\ and pi are however re- 
duced due to the removal of the singly-excited spin waves 
at the center of the ensemble. The concomitant reduction 
of the denominator of <? 2 (0) leads to an increase of the 
latter. For the longer ensemble (Fig. [3^b,d)), the center 
peak is not as pronounced in the simulation because of a 
higher average number of excitations, and is not evident 
beyond the statistical uncertainties of the data. 

Our numerical modeling employs Monte-Carlo simu- 
lations for 100 atoms randomly sampled according to a 
gaussian density with transverse waists w' x — w' y and 
longitudinal waist w' z , see Ref. [22j]. A single chan- 
nel model is used to calculate the interaction-induced 
detuning The annihilation fi,k and creation dj^ 

operators of the retrieved field determine g^ 2 \0) = 
(aj Co aj Co ak ak )/(aj Co ako) 2 - The retrieval process is de- 
scribed by the linear relation ak = -y/eSko + VI e £> 
where e is the overall efficiency for the atom-light map- 
ping, field transmission, and photoelectric detection. If 
the corresponding operator £ is assumed to be in the vac- 
uum state, the value of g^(0) is independent of linear 
losses and is equal to the corresponding spin-wave corre- 
lation function 5 f(0) = (StJtJ k< A )/(SlS k<) ) 2 . In 
the limit P « 1 the photoelectric detection probability 
P w e-q. 

The action of the probe beam is simulated by transfer- 
ring the atoms located in the illuminated part of the en- 
semble into a state that is not retrieved by the read pulse, 
and the spin-wave correlation function becomes spatially 
dependent. The calculated solid curves and shaded re- 
gions in Fig. 3 are in in good agreement with the ex- 
perimental data, which is achieved even though our sim- 
ulations do not account for the dynamical effects of the 
excitation blockade, suggesting that the dephasing of Ry- 
dberg spin waves is largely responsible for the observed 
features. 

In practice background light and detector dark counts 
contribute to the the photoelectric detection probabili- 
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FIG. 4: (color online), (a) <?' 2 '(0) as a function of photo- 
electric detection probability P that is varied by inserting 
partially transmitting optical filters into the retrieved field 
mode. The solid line is the function f + (sP(0) - 1)(1 - § ) 2 
derived by choosing the background mode described by an 
annihilation operator £ to be in a coherent state with average 
photon number B. Here P is a sum of contribution Pq due 
to the retrieved field and of background contribution B. (b) 
<7^(0) as a function of the probe beam position along the 
atomic sample z . Filled circles are for the case without sub- 
traction of background counts as in Fig. 3 and open circles 
are for the inferred values of </ 2 '(0) that would have been 
measured instead if there were no background counts. 



employed in Fig. 5 (a) and (b), respectively, for varying 
values of spin- wave storage time T s . Both g 2 (0) and ij 
decrease with T s reflecting the growth of the Rydberg 
interaction volume. The peak of g 2 (0) at z' « becomes 
prominent at T s ~ 0.8 /is when the spin- waves with more 
than two excitations have mostly dephased. 

We have demonstrated the emergence of localized col- 
lective atomic excitations in a strongly interacting atomic 
gas. In the future the extent of atomic correlations 
may be increased by using a smaller atomic sample and 
employing lasers with narrower linewidths. It is likely 
that the spin- wave coherence is considerably more sensi- 
tive than a single-atom coherence to various broadening 
mechanisms. Therefore imaging techniques of the type 
described in Refs. 3(3, 31 1 may be advantageous for ob- 
serving Rydberg crystals. 

We thank P. Goldbart, B. Kennedy, M. Pustilnik, and 
A. Zangwill for discussions. This work was supported by 
the Atomic Physics Program and the Quantum Memories 
MURI of the AFOSR and the NSF. 
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FIG. 5: (color online), (a) Intensity correlation function at 
zero delay, g' 2 '(0). (b) intensity of the signal, r\, both as 
a function of the probe beam position along the z-axis, for 
storage periods T a up to 0.8 [is, with w' x y = 8 fim, and w' z 
=42 fim. 



ties, so that the measured value g^(0) will differ from 
the quantity g s (0). To study the role of background 
photoelectric detection events in our experiment, we vary 
the probability of a photoelectric detection event P by 
inserting optical filters into the retrieved field with no 
probe beam used, Fig. 4(a). The data shows that over 
the dynamic range of our experiment the measured val- 
ues of g( 2 )(0) are not significantly affected by losses. To 
further assess the influence of the background counts, in 
Fig. 4(b) we display both g^(0) and the inferred quan- 

(21 

tity g s (0). The latter values are somewhat lower when 
the probe beam is near the center of the sample, with the 
same functional dependence. 

In the future it may be possible to observe the time 
evolution of the spin waves by achieving longer coherence 
times or a higher temporal resolution of the excitation- 
probe- retrieval sequence. As an example, we show g 2 (0) 
and rj normalized by its value when no probe beam is 
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